AD-A264  71 


'ATION  PAGE 


|  Pubiic  t< 
maintair 

suggestions  for  reducing  inisDu.ce!  i .  iu  «.c 
ana  to  me  Cfflce  of  Management  ana  Bud 

1  AGENCY  USE  ONLY  /Leave  blank) 


c-ur  per  response  rr^L.  j 
on  Send  com  merts  «•- ;  - 
,  i!efOr!ri?crn*,ationGix-»  *t 

*-aperwor»c  Reduction  Project  (Of 04  0168V  Warning’, on  [X~; 

2  REPORT  DATE 

March  1993 


4  TITLE  AND  SUBTITLE 

TERRAIN  MODELLING  USING  THE  SPLIT-STEP  PARABOLIC 
EQUATION  METHOD 

6  AUTHOR(S) 

A.  E.  Barrios 


7  PERFORMING  ORGANIZATION  NAME(S)  AND  AODRESS(ES| 

Naval  Command,  Control  and  Ocean  Surveillance  Center  (NCCOSC) 

RDT&E  Division 

San  Diego,  CA  92152-5001 


Professional  Pn|«T 

PR:  CDB0 
PE:  0602435 N 
WU:  DN48876U 


RMiNQ  ORGAN  VAI  i". 
T  NUMBER 


1 0  SPON  SORING  nnonroh  NO 
AGENCY  REPORT  NUMBER 


13  ABSTRACT  (Maximum  200  words) 


There  are  currently  several  software  programs  that  model  radiowave  propagation  over  irregular  terrain.  These  models 
use  a  combination  of  spherical  earth  diffraction,  multiple  knife-edge  diffraction,  wedge-defraction,  and  geometrical  optics  to 
arrive  at  a  solution  for  the  field  for  a  given  transmitter/receiver  geometry  and  a  specified  terrain  path.  More  recently,  para¬ 
bolic  equation  (PE)  methods  have  been  applied  to  model  propagation  over  terrain,  such  as  that  developed  by  Levy. 


9  3 


93-11408 


Published  in  Proceedings  IEEE  International  Conference  Radar  92 ,  No.  365,  Oct  1992,  pp  66-69. 


N$N  754001-2605500 


Stalwart  form  298  (BACK) 

UNCLASSIFIED 


\\,w  -•*>; ; ;  •  -v  •;•  v  v  *n\  •  .v  >  -  -v 


ICE  International  Conference  Radar  92 
12-11  October  1992,  Brighton,  England 


TERRAIN  MODELLING  US  INC  THE  SPLIT-STEP  PARABOLIC  EQUATION  METHOD 

Amalie  E.  Barrios 


Naval  Command,  Control  and  Ocean  Surveillance  center 
RDTiE  Division,  USA 


l-Accealon  Vor,  **’*!" 

^NTIS  CRA&J  & 

DTIC  TAB  A 

Unannounced  rj 

Justification 


INTRODUCTION 


By _ 

uistrfbution  / 

Availability  Codes 

Oist 

6d 

Avail  and/or 

Special 

There  are  currently  several  software  programs 
that  model  radiowave  propagation  over 
irregular  terrain.  These  models  use  a 
combination  of  spherical  earth  diffraction, 
multiple  knife-edge  diffraction,  wedge- 
diffraction,  and  geometrical  optics  to  arrive 
at  a  solution  for  the  field  for  a  given 
transmitter/receiver  geometry  and  a  specified 
terrain  path.  More  recently,  parabolic 
equation  (PE)  methods  have  been  applied  to 
model  propagation  over  terrain,  such  as  that 
developed  by  Levy  44-)". 

The  most  familiar  or  well-known  of  these 
models  is  the  Longley-Rice  model  (2).  This 
model  was  designed  for  low-altitude 
propagation  and  works  fairly  well  for 
diffraction  and  near-diffraction  regions.  A 
site-specific  propagation  model  for  general 
terrain,  called  SEKE,  was  devsloped  at 
Lincoln  Laboratory,  Ayasli  (3).  This  model 
is  based  on  the  assumption  that  the 
propagation  loss  over  any  path  (in  the 
frequency  range  from  VHF  to  X-band)  can  be 
approximated  by  one  of  the  multipath, 
multiple  knife-edge  diffraction,  or  apherlcal 
earth  diffraction  losaea  alone,  or  a  weighted 
average  of  these  three  basic  losses.  Another 
model,  developed  at  Ohio  State  University, 
Luebbers  (4),  is  based  on  the  geometricsl 
theory  of  diffraction  (GTD)  and  works  by 
determining  the  existing  rays,  for  a  given 
height/receiver  geometry  and  terrain  profile, 
from  a  family  of  16  ray  typea.  The  total 
field  at  the  target  is  then  found  by  adding 
the  ray  amplitudes  from  each  possible  ray. 

Each  of  these  models  has  various  limitations, 
but  the  main  limitation  they  all  share,  with 
the  exception  of  the  PE  model,  is  the 
inability  to  handle  ducting  or  non-standard 
range-dependent  environmental  conditions. 

SEKE  allows  •  variable  earth  radius  factor, 
but  this  assumes  a  constant  gradient  and 
horizontal  homogeneity.  Some  of  the 
literature  regarding  GTD  has  stated  that  this 
method  can  be  extended  to  inhomogeneous 
media,  but  this  author  has  not  seen  any 
published  reeults  for  such  cases. 


The  application  of  the  split-step  method  to 
model  tropospheric  radiowave  propagation  over 
the  ocean  has  been  well  documented  (6),  (7). 
The  problem  becomes  somewhat  more  complicated 
when  applying  this  same  technique  to  model 
radiowave  propagation  over  irregular  terrain. 

In  the  following  formulation,  the  atmosphere 
is  assumed  to  vary  in  range  and  height  only, 
makinq  the  field  equations  independent  of 
azimuth.  Beqinning  with  the  parebolic 
equation  derived  by  Fock  and  after  making  the 
envelope  transformation  (8),  the  equation 
that  must  ba  solved  is 

♦  2 ik.^t  *  kj(n'-l)+  -  0  (l) 

where  k ,  ia  the  free-space  wavenumber,  n  is 
the  index  of  refraction,  <t>  represents  a 
scalar  component  of  the  electric  field,  and  x 
and  z  are  the  spatial  cartesian  coordinates 
corresponding  to  range  and  height, 
respectively. 

For  propagation  over  terrain,  and  assuming 
horizontal  polarization,  equ.  (l)  is  subject 
to  the  range-dependent  boundary  condition, 
f(x,t-f(x))  -  0,  where  t(x)  is  a  general 
function  describing  the  terrain.  A 
transformation  ia  made  according  to  Beilis 
and  Tappert  (9),  which  generalizes  the 
"earth  curvature"  transformation,  and 
affectively  maps  the  range-dependent 
"terrain"  coordinate  system  to  a  flat  or 
smooth  earth  coordinate  system.  This  results 
in  a  "modified"  parabolic  equation  subject  to 
the  simpler  boundary  condition  that  the  field 
vanishes  at  the  surface  -  which  is  now  range- 
independent  in  the  new  coordinate  system. 

This  problem  can  then  be  easily  solved  by 
using  the  split-step  method  as  described  in 
(6). 

The  transformation  is  made  by  introducing  a 
change  of  variables.  Let 

X  •  x 

t  •  a-/(x) 


The  PE  modal  from  (1),  called  FDPEM  (Pinita 
Dlffaranca  Parabolic  Equation  Modal)  is  abla 
to  handle  range-dependent,  dueting 
conditions.  As  the  name  implies,  it  solve* 
the  parabolic  equation  by  using  finite 
difference  techniques.  However, 
compute tlOMlly,  finite  difference  methods 
can  ba  time  consuming, 

This  paper  presents  an  efficient  method  by 
which  one  can  determine  the  field  at  any 
point  above  the  earth's  aurfaca  in  tha 
praaeno*  of  ranga-dapandant  atmospheric 
conditions.  Tha  modal  is  basad  on  tha  split- 
step  Pourisr  algorithm  developed  by  Hardin 
and  Tapper t  (9)  te  solve  the  parabolic 
aquation.  Comparisons  are  mad*  against 
measured  data  and  tha  above  mentioned  models. 


where 

f(x)  •  t(x)-.£j 

and  defina  tha  scalar  component  of  the  field 
In  terms  of  the  new  coordinate  system: 

*(*,»)  ■  *(*,«)•'«•*.  <*> 

The  function  tfx)  describe*  the  actual 
terrain  and  can  ba  any  digitised  sat  of 
height/rang*  points.  xVi*  (whera  a  la  the 
earth's  radius)  takaa  into  account  the 
earth's  curvature. 

Substitution  of  equ.  (I)  into  equ.  (1)  yields 
the  modified  parabolic  aquation 


# 


1J  ■  (1 

,1?'  Tf.v 

«■’  -2g(t“  (x)  1  ) 

<1 


(3) 


*i  for  a  homogeneous,  •»  Inv.it.-.l  dm-i 
environment  over  Kast  Anql  la  .  I'rop.uj.ii  , 
loss  is  indicated  by  the  different  gray 
shades  and  the  antenna  is  1  orated  n>>t»r 

above  the  ground. 


with  the  new  boundary  condition,  \p(X,0)  «  0. 
t‘‘(x)  in  equ.  (3)  represents  the  second 
partial  derivative  with  respect  to  x.  Notice 
in  comparing  equ.  (3)  with  equ.  (1),  the 
inclusion  of  an  arbitrary  terrain  has 
effectively  produced  a  "new"  modified 
refractive  index.  This  is  consistent  with 
the  modified  refractivity,  or  M-unit, 
normally  used  in  tropospheric  wave 
propagation  over  the  ocean,  which  was  derived 
to  take  into  account  the  earth's  curvature. 
The  split-step  Fourier  method  is  then  applied 
to  give  the  solution  of  the  field  at  discrete 
range  steps  for  all  target  heights  under 
consideration. 


Since  equ.  (3)  does  not  allow  for  prop.iqa*  i  *s. 
over  vertical  obstacles,  such  as  cliffs  oi 
buildings  (  t''(x)  will  be  undefined  for  smti 
cases  ) ,  a  special  case  can  be  made  by 
simply  eliminating  the  field  immediately 
adjacent  to  such  obstacles  and  propagating 
the  field  forward  as  usual.  This  docs  not 
violate  any  conditions  in  the  split-step 
model  as  the  PE  approximation  inherently 
neglects  backscatter.  Figure  6  shows  a 
coverage  diagram  for  such  a  case.  Excellent 
agreement  was  found  when  compared  against 
FDPEM. 

CONCLUSIONS 


This  method  offers  a  numerically  efficient, 
full  wave  solution  to  the  field  because  of 
the  implementation  of  the  Fast  Fourier 
Transform  (FFT)  in  the  computer  model. 
Computer  execution  times  increase  for 
increasing  frequency  and/or  large  propagation 
angles,  i.e.,  steep  terrain  slopes.  However, 
the  split-step  method  remains  more  efficient 
than  finite  difference  techniques  for  these 
extreme  cases. 

RESULTS 

Figures  1  and  2  show  comparisons  between 
measured  data  and  the  Terrain  Parabolic 
Equation  Model  (TPEM) .  The  results  are 
displayed  as  height  versus  propagation  factor 
(or  field  strength  relative  to  free  space)  in 
dB,  with  "height"  referring  to  height  above 
the  ground  at  the  particular  receiver  range 
shown.  The  terrain  profile  is  also  shown  for 
each  case.  For  both  figures  the  environment 
was  assumed  to  be  a  homogeneous  standard 
atmosphere.  The  data  shown  in  Fig.  1  were 
taken  from  ref.  (2).  Results  from  the  model 
SEKE  are  also  shown.  Data  from  ref.  (3)  are 
shown  in  Fig.  2  along  with  results  from  the 
GTD  and  the  Longley-Rice  model.  In  both 
figures  TPEM  agrees  very  well  with  the 
measured  data  and  the  other  models. 

In  1946,  the  Navy  Electronics  Laboratory  (now 
Naval  Command,  Control  and  Ocean  Surveillance 
Center)  conducted  an  experiment  over  the 
Arizona  desert  to  study  atmospheric 
inhomogenelty  and  the  irregularity  of  the 
terrain  (10).  The  transmitter  was  located  at 
Gila  Bend  with  a  receiving  station  at 
Sentinel.  This  site  was  chosen  because  of 
the  large  diurnal  variation  in  surfaca 
temperature  common  in  the  daeert,  leading  to 
strong  inhomogeneous  ducting  conditions. 
Meteorological  measurements  were  teken  at  the 
transmitting  and  receiving  station!,  and  at 
two  stations  along  the  path.  The  terrain 
profile  end  the  meteorological  measurements 
for  0300  February  6,  1946  (plotted  as  height 
vereua  M-unite)  ara  ahown  in  Figura  3. 

Arrowa  Indicate  the  location  of  the  measuring 
atatlona  and  tha  "0"  haighta  on  tha 
rafractlvity  acalas  correspond  to  tha  local 
ground  aurfaca.  Figura  4  shows  the  measured 
data  along  with  that  predicted  by  TPEM  for 
both  tha  aetuai  terrain  and  tha  smooth  earth 
case.  Aa  tha  figure  shows,  one  must  also 
taka  into  account  tha  aetuai  tarraln  aa  wall 
aa  tha  environment  to  make  proper  field 
pradiotlons. 

To  demonstrate  the  power  of  the  epllt-etep 
method,  a  coverage  diagram  la  shown  in  Figura 


A  numerically  efficient  method  has  been 
■presented  to  model  tropospheric  radiowave 
propagation  over  irregular  terrain  in  the 
presence  of  range-dependent  non-standard 
environmental  conditions.  Results  from  this 
model  were  compared  against  measured  data  and 
other  existing  models  and  was  shown  to  give 
excellent  agreement.  This  work  is  in  the 
public  domain. 
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Figure  1.  Comparison  between  TPEM,  SEXE,  and 
measured  data  at  167  MHz  for  standard 
atmosphere  over  the  terrain  shown. 
Transmitter  height  is  at  18.3  meters  above  the 
ground  at  0  range.  Receiver  range  is  at  35.5 
km. 


Figure  2.  Comparison  between  TPEM,  GTD, 
Longley-Rice,  and  measured  data  at  110.6  MHz 
for  standard  atmosphere  over  the  terrain 
shown.  Transmitter  height  is  at  1.36  meters 
above  ground  at  0  range.  Receiver  range  is  at 
6.6  km. 
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Figure  3.  Terrain  and  environmental  profiles 
for  Figure  4.  Arrows  indicate  location  of 
meteorological  measurements. 


Range-dependent  Ducting  Atmosphere  -  Gila  Bend  to  Sentinel 
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Figure  4.  Comparison  between  TPEM  for  actual 
terrain,  TPEM  for  smooth  earth,  and  measured 
data. 


